Solar system constraints on asymptotically flat IR modified Hofava gravity through 

light deflection 

Molin Liu 1 !*! Junwang Lu 1 , Benhai Yu 1 , and Jianbo Lu£_| 

1 College of Physics and Electronic Engineering, Xinyang Normal University, Xinyang, 464000, P. R. China 
2 Department of Physics, Liaoning Normal University, Dalian, 116029, P. R. China 

In this paper, we study the motion of photons around a Kehagias-Sfetsos (KS) black hole and 
obtain constraints on IR modified Hofava gravity without cosmological constant (~ Aw)- An 
analytic formula for the light deflection angle is obtained. For a propagating photon, the deflection 
angle 8ip increases with large values of the Horava gravity parameter u. Under the UV limit 
oj — > oo, deflection angle reduces to the result of usual Schwarzschild case, 4GM/R. It is also 
found that with increasing scale of astronomical observation system the Hofava-Lifshitz gravity 
should satisfy \ojM 2 \ > I.1725x 10" 16 with 12% precision for Earth system, \ojM 2 \ > 8.27649 x 10~ 17 
with 17% precision for Jupiter system and |o;M 2 | > 8.27650 x 10~ 15 with 0.17% precision for solar 
system. 
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I. INTRODUCTION 



It is known that the Lagrangian of original Lifshitz scalar field theory is 

C -- 



d 2 xdt 



(1) 



which has a line of fixed points parameterized by Xl with anisotropic scale invariance ll]. So such fixed points 
are the so-called Lifshitz points. Recently, Hofava adopted one Lifshitz scale in time and space, and obtained a 
renormalizable gravity theory at these Lifshitz point [2H4|. It is known also as Hofava-Lifshitz (HL) gravity. It 
exhibits a broken Lorentz symmetry at the short distances, but it reduces to usual General Relativity (GR) at the 
large distances due to these non-contribution of higher derivative terms with particular A = 1. The parameter A 
controls the contribution of the extrinsic curvature trace. However, it is possible that HL at large distance could 
deviate from GR for the incomplete diffcomorphism invariance. In this paper, we only consider the simple case of 
A = 1. Since HL gravity theory _was put forth, HL gravity is intensively investigated in many aspects. They are 



developments of basic formalism [5 - |24l. cosmology 



solutions and their thermodynamics 



37 



25 



dark energy and dark matter 



37 



various black hole 



- 5r3 | and so on. 

With the subsequent developments of the HL gravity, the following problem raised in front of us: how does the matter 
fields influence the geometry? Several researchers made efforts to find the analogue of the matter energy-momentum 



tensor acting gravitational source as GR. The study of its geodesic in HL theory is the first step for finding a solution. 
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Hence, some pioneering works addressed this aspect with various methods. They are the optical limit of a scalar held 



theory 



65], super Hamiltonian formalism foliation preserving diffeomorphisms [67[, Lorentz- violating Modified 



GO 



62 



64]. In particular, in 



Dispersion Relations [68| as well as observational kinematic constraints below mentioned 
the framework of the optical limit Capasso and Polychronakos [65J gave us a comprehensive study of particles in the 
HL gravity. Based on the optical limit, the deformed geodesic equation was derived from a generalized Klein-Gordon 
action. Such results clearly show that the fact that the particles move along geodesic, which depend on the mass of 
particles, is not granted. Deviations from geodesic motion appear both in flat and Schwarzschild-like spacetimes. This 
type deformed kinematics could allow the presence of superluminal particle for particular A ( or A's in Ref. 65|) which 
could not give any restriction on the velocity. The particles with given A's could be always accelerated to any velocity 
under an action of constant force. Else, the higher is the energy of the particle, the larger will be the deviation from 
the geodesic motion determined by the low-energy four dimensional metric g^ v . Similar results deviations from GR 



also are found in 



with various above-mentioned approaches. 



On the other hand, a specific case of asymptotically flat IR modified HL gravity under the cosmological constant 
limit Aw — > in HL gravity is proposed by Kehagias and Sfetsos [53 1 . It reduces to GR gravity in IR regime according 
to coupling A. The massless excitation with two propagating polarizations emerges in the Minkowski vacuum, which 
is admitted naturally, with the transverse traceless part of the perturbation. At the most important, one static 
spherically symmetric black hole solution of HL gravity i.e. KS black hole is obtained, which can be treated as 
counterpart of the Schwarzschild case in GR. Comparing its geometry with GR, it is also found that this solution 
has the same Newtonian and post-Newtonian limits which ensures its soundness when it is confronted with the 
astronomical observational data. Since then, there are some works involving the phenomen ology o f the deformed HL 
black hole, such as strong field gravitational lensing 0,[58|, black hole's quasi-normal modes 58, 5tj], timelike Geodesic 



Motion 



601 ] . thin accretion disk [611 ] . In Konoplya's work 58[, the gravitational lensing and quasi-normal modes are 



considered, which are potentially observable properties of black holes beyond the post-Newtonian corrections. It is 
found that the bending angle is smaller and the real oscillation frequency of quasinormal modes (QNM) is larger than 
that of GR case. The QNM can be lived longer in this deformed HL black hole. In Chen's work [57j, the angular 
position and magnification of the relativistic images are obtained with strong field gravitational lensing. Unlike the 
former Konoplya's work [5^ , a scheme to distinguish a deformed HL black hole from a Reissner-Norstrom black hole 
is introduced. Comparing with the latter work, the angular position decreases more slowly with bigger HL gravity 



parameter u). Else, angular separation increases more rapidly. In Chen's another work [591 ] . it is also found that the 
scalar perturbations decay more slowly in the HL gravity. In Ref. [soj l , th e time-like geodesic motion is analyzed with 



the effective potential in this deformed HL gravity. In Harko's work [61|, by using the accretion disk properties, one 
can obtain some important properties such as the energy flux, the temperature distribution, the emission spectrum 
and so on. Unlike the usual GR, there are some particular signatures appeared in the electromagnetic spectrum which 
maybe supply a observable path for HL gravity. 

With the development of KS black hole, there are some works referring to field parameters constraints such as the 
solar system orbital motions 62], the extrasolar planets and double pulsar 63(, the solar system three classical tests 



64| and so on. In 62|, |63[ lower bounds on uj are obtained via the post-Newtonian motion equation. 



As we know that, the asymptotically flat HL gravity with A = 1 reduces to the Schwarzschild space with the limit 
of cj — > +00. So in the UV limit, the field parameter u should be large enough to ensure the UV Horava gravity is in 
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conformity with GR gravity well. Motivated by above various factors, one aim of this paper is to find the lower limit 
of HL gravity field parameter u. Considering the last work [64j |. the numerical constraints results are obtained directly 
from integrating metric function by using the static and spherically symmetric metric. Hence, another aim of this 
paper is to find the analytic expression of light deflection. We adopt one approach of Lagrange analysis and study the 
light deflection in the KS black hole in HL gravity under the weak field and the slow motion approximation. Then, 
we also compare our results with the observational data (long baseline radio interferometry, Jupiter measurement, 
Hipparcos satellite) and we also obtain corresponding constraints on HL field parameter oj for solar system, Jovian 
planet system and Earth system, respectively. 

This paper is organized as follows. In section II, we present the action of HL gravity and spherically symmetric 
Kehagias and Sfctsos black hole. In section III, we calculate the general solution to motion equation of photons in 
deformed HL gravity. In section IV, we use astronomical observations to constrain the HL gravity. In Section V we 
compare our results with other constraints obtained from inside and outside solar system. Section VI is devoted to 
the conclusions. We adopt the signature (— , +, +, +) and put %, c, and G equal to unity. 



II. THE ASYMPTOTICALLY FLAT BLACK HOLE SOLUTION IN THE DEFORMED HL GRAVITY 

In this section, we review briefly the KS black hole solutions under the limit of A^y — > with running constant 
A = 1 in the IR critical point z = 1. The spacetime geometry is, in the ADM formalism \&\ . 

ds 2 = -N 2 dt 2 + g l} (dx l + Wdt) (dx j + N 3 dt) , (2) 

where the classical scaling dimensions of the field are 

[ffy]. = 0. [N}s=0, [N% = 2. (3) 

The action for the fields of HL theory is 



S = J dtd^gN^ (K ijK - XK 2 ) - £-(■„(■>■> + ge^V.i^ 

-^f^R^ + (i? (3) ) 2 + A^ (3) - 3A^) + /i^)}, (4) 

where the second fundamental form, extrinsic curvature K^j, and the Cotton tensor C lJ ' are given as follows, 

C ij = e rkl^j k { R <J!.)j _ IflW^ . (6) 

Here, k, A, w are dimensionless coupling constants and fj,, Ayy are dimensional parameters having dimensions of 
masses [//] = !, [A w ] = 2. The last term of metric (jl]) represents a soft violation of the detailed balance condition. In 
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the limit of A\y — ► 0, we can obtain a deformed action as follows, 

S = J dtd^xiCo+d), (7) 



Co = ^gN^iKijK* -\K 2 )}, 



(8) 



A - V9N{ ~*S 2 ~ £ ( ~ f - fgfl* ) + (9) 



32(1 - 3A) 2w 4 V 2 / V 2 

Comparing these with the usual GR ADM formalism, we can obtain the speed of light c, the Newton's constant G 
and the cosmological constant A 



For the particular case of A = 1 with u> = 16/j, 2 /k 2 , a spherically symmetric black hole solution is presented by 
Kehagias and Sfetsos 53j, which also is corresponding to an asymptotically flat space, 



ds 2 = -f(r)dt 2 + ^r+r 2 (d6 2 + sin 2 6d<p 2 ) . (11) 
J{r) 



The lapse function is 



f(r) = 1 + tor 2 - y/r (uV + 4wM), (12) 



where the parameter M is an integration constant with dimension [M] = — 1. Using the null hypersurface condition, 
one can find there are two horizons, inner r_ and outer event horizon r + in this space, 



2uM 2 

The Hawking temperature is 



r± = M[l±dl--^ . (13) 



T, = (14) 

The corresponding black hole thermodynamics and its generalized uncertainty principle are studied in Ref. [jf]]. It 
is also found that if r 3> (M/w) 1 / 3 , this KS solutions will reduce to usual Schwarzschild case and the lapse function 
can be rewritten as 

/«!- — + 0(r" 4 ). (15) 
r 

Apparently, the former integration constant M is the real mass of black hole. 



III. MOTION OF PHOTONS IN HL GRAVITY: GENERAL SOLUTION 

By using the lapse function of KS black hole (TTT|) , the usual Lagrangian of test particle is given as 

2C P = -f(r)i 2 + j^ +r 2 (e 2 + sin 2 ^ 2 ) , (16) 
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where the overdot represents differentiation with respect to the afHne parameter £ along the geodesies. The Euler- 
Lagrange equations below give us the kinetic equation of photons, 

Without loss of generality, it is assumed that the observer is confined to orbits with 9 = tt/2 and 8 = 0. Hence, the 
Lagrangian C p reduce to 

2£ p = -f(r)i 2 + ^+r 2 ^ 2 . (18) 

According to the Euler-Lagrange equations (|17p and reduced Lagrangian C p (JTHJ) , we can identify two constant motions 
for cyclic coordinates ip and t, 

f(r)i = E, (19) 
r 2 p = L, (20) 

where E and L are the linear momentum and angular momentum at infinite distance. According to the condition for 
photons g^vP^p" = where p M = = is the photons momentum, we can obtain a relation of parameters (i, f, ip) 
as, 

-f{r)i 2 + ^-t 2 + r 2 p 2 = 0. (21) 
Dividing Eq. (|2ip by ip 2 , we can obtain motion equation as 



dt \ 2 1 / dr x 2 



Substituting Eqs. (jT9|) . (1201) into Eq. (|221) . we can obtain the final photons' orbital equation 

(23) 



1 dr\ 2 1 1 



j 2 dip J b 2 B 2 {r) 1 

where E and L are coupled via 6 = L/E which represents the effective sighting range or impact parameter. The 
corresponding photon effective potential 1 / B 2 (r) is given as 

B 2 (r) = /" =. (24) 

1 + ujr 2 - \Jr (u 2 r 3 + AujM) 

Fig- (H} clearly illustrates that the potential becomes smaller for a larger oj. One also could find that with the bigger 
HL gravity parameter ui the curves tends toward the usual GR Schwarzschild result which is drew with bigger red 
dotted line in Fig.([T]). This phenomenon is in agreement with the Schwarzschild limit r 3 ^> M/u) in which for 
the fixed variables r and the central mass parameter M, the rewritten Schwarzschild limit is w > M/r 3 . Hence, the 
bigger lu will lead to a final Schwarzschild case. 

The constant impact parameter b in Eq. (|23p can be eliminated via calculating the derivation with respect to (p. 
So the photon Binet equation can be obtained 

^ =-u + 3u 2 V^M 2 ~(loM 2 +4u 3 V\ (25) 
dp 1 ' 
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0.04 




FIG. 1: The effective potentials 1/B 2 (r) of photons with unit mass and various parameters lo. The usual Schwarzschild case 
in GR is denoted with bigger dotted red line. 



where u — M jr. Considering the fact of that all astronomical observations can be explained very well by the 
Schwarzschild gravity in solar system, we adopt the Schwarzschild approximation r 3> (M/uj) 1 ^ 3 . Hence, the final 
general solution to photon orbital equation (|23|) is written directly as 



u = uq cos ip + Uq (l + sin 2 p) + — ( tp simp — — cos 3p — cos 5ip ] . (26) 

2 \ 8 120 J 

Here, the parameters uq and a are uo = GM/R and a = — 15uq/4ujM 2 where R is the solar radius and M is the 
solar mass. In the right hand side of Eq. (|26p . the first term u$cosp is a straight line normal to pole axis without 
considering the influence of gravity, the second term w 2 , (l + sin 2 ip) is the standard form coming from general relativity 
gravity, and the last one ^ (ipsiiKp — i cos3p — cos5y) is an improvement of HL Gravity. Part of some related 
calculations are given in the following section. 



IV. SOLUTION TO PHOTON ORBITAL EQUATION WITH HOrAVA-LIFSHITZ GRAVITY 

Under Schwarzschild approximation r ^> (M/w) 1 / 3 , the photon Binet equation Eq. (l25T) reduces to a solvable form, 

d ^-- U + 3U2 ~ 6 ^ (27) 

Considering the small u, the high-order small quantities u 2 and u 5 can be ignored firstly. So we can obtain the 
zcroth-order approximate solution 

u = u Q cos tp, (28) 

where parameter uq — GM/R is a constant about solar system. Mathematically, this solution is a straight line 
normal to the pole axis. In order to simplify the calculation, we only replace the high-order terms with the zeroth- 
order solution in Eq. (|2"T|) not the fussy iteration step by step. Hence, one can get 

—j = -u + 3uq cos 2 p -pr cos 5 p. (29) 

dp* ujM^ 
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(30) 



After decreasing the powers of cos 5 ip term, Eq. ([29|) can be reduced a new form 

d 2 u 
dip 2 

Because the right hand of Eq. (j3"0)) includes four terms, we can decompose this ordinary differential equation (f3"U|) into 
corresponding four sub-equations to solve. Here, we list these particular solutions in the Table [J] 



o 2 2 3u o r 15u o o 15w o 

3U ° C0S ^-^ COS5 ^-8^ COS3 ^-4^ COS ^ 



TABLE I: The particular solutions to Eq.pO 



Sub-equation 



Particular solution 



^ + u = 3u 2 cos 2 ip 



u — Uq (l + sin 2 ip) 



d~ a 
d s - 



jk cos 5(p 



2io cos5 V 



u = — ^ COS 3^2 



a cosv? 



Hence, the total particular solution is 



u 



4(1 



, s a / 1 1 

+ 2 I ^sin<^ — -cos3<^- Y^7j cos5< ^ 



And the first-order approximate solution to Eq. (|27p is given as 

a 



U = Uq COS <p + Uq (l 



sin 2 



y> sin y> - 



— cos 3<p cos hip 



(31) 



(32) 



2 V 8 120 

The azimuth angles of zeroth-order approximate solutions (|28j) are ±7r/2 at very great distance. However, the 
azimuth angles of the first-order ones (1321) are ±(7r/2 + /3) in u = 0. The deflection angle /3 is a small quantity which 
satisfies 



u sin/3 + Uq(1 + cos 2 0) + - 



| + cos ,3 - - sin 3/3 + — sin 5/3 



0. 



(33) 



We use the Taylor expansions of sine and cosine functions and keep the basic term. The deflection angle is written 
formally as 



Sep = 2/3 



R V RM + 8RM J ' uj + \ ARM 3 + 32RM 3 J ' w 2 + [Zp 



(34) 



where the first term AGM/R is the usual result of the pure Schwarzschild spacetime |7l|. The other high-order 
correction items containing parameter oj are the contribution of HL Gravity. It affects the light deflection only 
through the particular solution (|31[) . but the zeroth order approximate (|28[) is unchanged. Furthermore, the well 
asymptotically behavior of deflection angle, Sip — > AGM /R under the limit cj — > +00, could ensure UV HL gravity 
does not break the classical GR gravity. The deflection angles versus 10M 2 are illustrated in Fig|2]in which Sip increases 
with bigger uM 2 . Furthermore, the upper limit of this curve is the usual Schwarzschild case i.e. AGM/R. 



CONSTRAINTS ON HOrAVA-LIFSHITZ GRAVITY FROM ASTRONOMICAL OBSERVATIONS 



In this section, we will use the observational data to confine the Horava-Lifshitz Gravity based on the deflection 
angle Eq. (|3"4"f . It should be noticed that in order to analyse exactly the results Eq. (|3"4l , we move from the natural units 
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4GM. Ft 




-2 

1(T ! 10~ 1 10° 10 1 10 2 10 s 

toM 2 



FIG. 2: Deflection angle Sip versus HL Gravity parameter loM 2 with uo = 1- The solid bold line denoted by AGM/R is the 
usual Schwarzschild case [57]. 



used so far to the SI system by restoring c and G throughout the relevant formulas. In order to obtain experimental 
constraints from the light deflection result, the final light deflection angle 5(f should be expressed in terms of the 
deviation A^rj from the general relativity prediction StpaR for the Sun. 

5<p = 6(p GR (l + A LD ), (35) 

where SipaR — AGM / c 2 R. The best available constraints on SipaR come from long-baseline radio interferometry which 



shows that |A^n| < 0.0017 by Robertson et al. [72j and Lebach et al. [73[. Submitting the deflection angles Eq. (j3"4")l 
into Eq.p5p and neglecting the higher-order terms, we can therefore infer an lower limit of Horava-Lifshitz gravity 
parameter 

|wM 2 | > 8.27650 x 10~ 15 , (36) 

for the solar system. However, it is important to bear in minds that the parameter uj characterizing the HL gravity 
metric (jlip is not a naturally universal constant like G or c, but it may in principle vary from one space to another. 
HL gravity as an alternative to 4D general relativity is therefore best constrained by the application of two or more 
tests to the same system. 

With this in mind, according to the deviation |Ald| < 0.17 from the Jupiter measurement of light deflection by 
Treuhaft and Lowe [7J|, we can obtain 

\ujM 2 \ > 8.27649 x 10~ 17 . (37) 

It has also been proposed by Gould [75I to measure the light deflection due to the Earth by using the Hipparcos 
satellite, with an estimated precision of 12%, which leads a sensitive values for the Earth system 

\ujM 2 \ > 1.1725 x 10" 16 . (38) 
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VI. COMPARISON WITH OTHERS CONSTRAINS FROM SOLAR SYSTEM AND EXTRASOLAR 

PLANETS 



Here we compare our results with the constraints obtained by other researchers using solar system and exoplanetary 
data. The solar system tests include perihelion precessions 



G2. 64], rada.r echo delays and light bending 



6J. The 



cxtrasolar planets have been recently used bylorio and Ruggiero [63j as well. 



64J , the constrains are obtained via the analysis of arbitrary spherically 



In the former classical solar system tests 
symmetric spacetimes [64| . which is different from the dynamic Lagrangian analysing in this paper. It should be noted 
that the parameters uiq and ipo used in 62H64} coincide with our ujM 2 . For the perihelion precession case, the circular 



orbit uq is determined by the nonlinear algebraic equation (51) in Ref. [64J which could bring us numerical solutions. 
So using the planet Mercury data, one can obtain the corresponding result ujM 2 ~ 7 x 1CP 16 with the 43 arcsec per 
century of precession angle immediately, which is in the same order with our result of Hipparcos-Earth system, i.e. 
1.2 x 1CP 16 . For the radar echo delays case, the time interval between emission and return receiver is determined by 
the integration of lapse function, 

At = J h | e [MVra ? )-(NA i ^)]/2_ 1 |. (39) 



Though this general expression is simple, it is hard to solve and obtain a explicit form of time delays, except for the 
numerical solution, since there is a square root term in the KS lapse function. The numerical integration shows the 
result luM 2 ~4x 10 -15 which has an exact identical order of magnitude with our long-baseline radio interferometry 
result, i.e. 8.27650 x 10~ 15 under analytic analysis perspective. For the light deflection case, it is slightly emphasized 
that, in spite of considering light bending both in Ref. 64J and this paper, the based methods have nothing in common. 



Furthermore, the analysis of light bending based on the clear analytic expression of deflection angle (|34|) is grasped 
in this paper. Then, we have a look at the relating results in Ref. [64] . The photon's motion equation satisfies the 
integral expression (55) of Ref. [64j . Hence, the deflection angle is Acf> = 2\</)(xq) — <f>(oc)\ — it. Like the former radar 
echo delays case, since there are some square root terms, the numerical solution is obtained rather than analytic one. 
Submitting the classical observational value A<f> = 1.7275 arcsec, the result shows parameter luq — 10 -15 , which is 
the same order as our results ujM 2 ~ 8.3 x 10 -15 based on the analytic perspective. Meanwhile, this results of the 
constraints on asymptotically flat IR modified Horava gravity are compatible with the constraints in Jupiter and Earth 
systems. The corresponding compatibilities of parameter uiM 2 are in the orders of 1 x 10~ 15 (Harko and Kovacs) > 
8.27650 x lCT 15 (Solar system) > 1.1725 x 10~ 16 (Earth system) > 8.27649 x 10~ 17 (Jupiter system). Obviously, the 
constraints on Horava gravity field parameter lo are more restrictive via the analytic method of this paper. 

Then we consider the orbital motions tests 62, 0] in the following words. It should be noticed that in the 
approximation of weak field and slow motion, orbital motion mainly refers to the secular precession of the longitude 
for the pericenter test particle which has not perihelion at all. So orbital motion is different from that of the general 
rclativistic Einstein's pericentre precession. The constraints of orbital motion are divided into two parts, inside 
and outside 63j our solar system. The secular precessions of the longitude of the pericentre w are calculated for a 



test particle by using the Gauss and Lagrange perturbative approaches. Then, based on the EPM2008 ephemerides 
constructed from a data set of 550000 observations ranging from 1913 to 2008 [t3 | . they compared their results with 
the corrections Azu to the standard Newtonian/ Einsteinian planetary perihelion precessions and found the lower 
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bounds on KS solutions. For the case of solar system planets [62j, inferior planets (Mercury, Venus, Earth, Mars) and 
outer planets (Jupiter, Saturn, Uranus, Neptune, Pluto), the lower limit of uiM 2 is 



3(GM) y / 2 (4 + e 2 ) 



2[Ari7±5(Ari7)] MaxC 6a9/2(l- e 2)3 



(40) 



where a is the semi- major axis and e is the planet's eccentricity. So the corresponding constraints on HL gravity 
ranging from | i/>q | > 10~ 12 to \i/jq\ > 10 -24 are shown in following tables |62j. Interestingly, the larger value 8 x 10 -12 



in their work |62| coming from Mercury is about 1000 times more than that of our result Eq. p6[) coming from radio 
intcrferometry. For the same Earth, we also find that Iorio and Ruggiero's result 6 x 10 -14 is almost 513 times more 
than that of our result Eg. (13"B"1) . In contrast, our result Eq. (j3"7} is almost 30 times more than Iorio and Ruggiero's 
result 3 x 10 -18 for the same Jupiter plant. 



TABLE II: Lower bounds on |^o| of of the inner planets and outer planets in solar system 



inner 


Mercury 


Venus 


Earth 


Mars 




planets 


8 x 10" 12 


6 x 10" 14 


6 x 10" 14 


8 x 10" 15 




outer 


Jupiter 


Saturn 


Uranus 


Neptune 


Pluto 


planets 


3 x 10" 18 


1 x 10" 18 


2 x 10" 22 


2 x 10" 23 


8 x 10" 24 



Then we consider the other orbital motions type constraints coming from the extrasolar planets shown also by 
Iorio and Ruggiero [63| . Like their early solar system work [62j , the same standard Gauss perturbative approach was 
adopted. The correction to the Keplerian period is 

_ ^{GMf/ 2 {l + \e 2 ) 
u ° Wo c 6 a 3 / 2 (l-e 2 ) 5 / 2 ' [ ' 

where ujq = uM 2 = i/jq. By using the limit of e — > 0, one can equates the centripetal acceleration to the Newton + 
KS gravitational acceleration for a circular orbit. The final expression of ujq in HL gravity is 

4tt(GM) 5 / 2 

= -£»&r, (42) 

where d represents the current fixed radius of the circular orbit and AP is the different value of correction orbital 
period to the third Kepler law. Comparing these with the phenomenologically determined orbital periods of the 
transiting extrasolar planet HD209458b "Osiris", the result uj > 1.4 x 10 -18 which is about one fiftieth of our 
minimum Jupiter planet value Eq. (|37p . Obviously, for the extrasolar planets give us a more rigorous constraints. 



VII. CONCLUSION 



In this paper, we have studied the constraints on HL gravity from light deflection astronomical observations including 



72 



73 1 , Jupiter measurement 74] , Hipparcos satellite [75] . Here we summarize our 



long-baseline radio interferometry 
findings. 

We need briefly specify the parameters A and Aw firstly. A is dimensionless coupling constant which controls the 
contribution of the extrinsic curvature trace. In this paper, we only consider the specific case of A = 1 which is 
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corresponding to usual GR case at large distances. The parameter Aw is the cosmological constant of HL gravity. 
For the various HL gravity models with Lifshitz points, Aw has the different reduced values. For example, for Z = 1 
HL gravity in this paper, its value is Aw = 3 A/2 where A is the cosmological constant of usual 4D GR. For Z = 4 
HL gravity, its values are Aw = 3 A/2 for (3 + 1) dimensions and Aw = A for (4 + 1) dimensions Q. Considering 
the asymptotically flat IR modified HL gravity, we adopt Aw — in this paper. 

By using the conventional Lagrange approach for the calculus of variations, the the motion of photons is studied 
in the deformed HL gravity. The photons' orbital equation (|23|) is coupled with the key HL Gravity parameter uj in 
the potential function 1/B 2 (r). Hence, a scattering problem in a central force field occurs here. It is found that the 
photon collisional potential, which is larger for smaller u) is larger than in the usual GR case. Furthermore, for large 
values of uj, the HL gravity reduces to the Schwarzschild case which is also justified by the Kehagias and Sfetsos's 



Schwarzschild limit 



53|. 



All the classical tests in the solar system are in agreement with GR, especially those by long-baseline radio inter- 
Jupiter measurement [74] , Hipparcos satellite [isl ■ So in this paper we use the Schwarzschild limit 



ferometry 



72 



53] and the difficult Eq. (|25p is directly solved according to the linear superposition principle of ordinary differential 



equation. In the final result of Eq. (|34[) , wc find that the deflection angle value increases monotonically for large values 
of the HL gravity parameter ujM 2 (or uj). For a small magnitude of uj (about 10~ 17 ~ 10~ 15 here), the modification 
of HL gravity is dominated by the second term (~ 1/oj) according to Eq. (p4)l . So the result of HL gravity is less than 
the usual GR value AGM/R. This relation of HL and GR is also clearly illustrated in FIGEJ The similar behavior 
will be found in the strong field gravitational lensing 571 . 1581 . Finally, the constraints on HL gravity from long- 
baseline radio interferometry 72, 73 1, Jupiter measurement [zj], Hipparcos satellite 3] are |wAf 2 | > 8.27650 x 10~ 15 , 
|wM 2 | > 8.27649 x 10~ 17 and \ujM 2 | > 1.1725 x 10~ 16 for Solar system, Jupiter system and Earth system, respectively. 

For the comprehensive consideration of the constraints on asymptotically flat IR modified HL gravity, we also 
compare our results with others similar works including usual GR tests [64|, orbital motions inside 



62J and outsider 



63j solar system. Firstly, for the comparison with usual GR tests 64|, we find our constraint from Hipparcos-Earth 
system (~ 10~ 16 ) is concordant with the result of perihelion precession of Ref.jsjj]. Meanwhile, our constraint from 
long-baseline interferometry (~ 10~ 15 ) is in agreement with their result of radar echo delays case. However, for the 
same light deflection test, our results could give a more lower bound than their value (~ 10~ 15 ). Secondly, for the 
comparison with orbital motions inside solar system [62I ]. we find that our result (~ 10~ 16 ) gives a lower bound, 
about two orders of magnitude, than that of their result (~ 1CP 14 ) for the same Earth system. However, our result 
(~ 1CP 15 ) gives a higher bound, about three orders of magnitude, than that of their result (~ 10~ 18 ) for the same 
Jupiter system. Thirdly, for the comparison with orbital motions outsider solar system the lowest bound from 
our Jupiter case (~ 1CP 17 ) is higher one orders than that of HD209458b "Osiris" (~ 10~ 18 ). In general, according to 
the above situation, the HL gravity field parameter ujq or tpo should be in the range of [10 -12 , 10 -24 ]. The constraint 
from solar system mainly concentrates in the range of [10 -17 , 10 -15 ], except for the orbital motions. On the contrary, 
the constraint from orbital motions distributes very separately, in which the largest bound (~ 1CP 12 ) and the lowest 



12 

bound (~ 1CP 24 ) come from in Mercury system and Pluto system, respectively. 
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